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Abstract 

Solution of the nonlinear Klein-Gordon equation perturbed by small 
external force is investigated. The perturbation is represented by finite 
collections of harmonics. The frequencies of the perturbation vary slowly 
and pass through the resonant values consecutively. The resonances lead 
to the sequence of the wave packets with the different fast oscillated car¬ 
riers. Full asymptotic description of this process is presented. 

Introduction 

In this work we study the problem of a generation of sequences of solitary 
packets with different carriers in the optical fiber. The nonlinear Klein-Gordon 
equation is studied as a modeling equation. We consider this equation which is 
perturbed by a small external driving force with finite collection of modes. The 
wave packets appear due to passage through a resonance by different modes of 
the external force. After the passage through the whole resonances the solution 
contains the full collection of the solitary packages of waves with different car¬ 
riers. The envelope functions of these packages satisfy to nonlinear Schrodinger 
equation (NLSE). 

In general the derivation of NLSE for small solutions of nonlinear equations 
is well known p m m and justified [J. Our solution has a more complicate 
structure. Before all resonances the solution has an order of the perturbation 
and defined by the external force. After the passage through the resonance of 
the last mode of the perturbation the solution has the order of the square root 
of the order of the perturbation and satisfies NLSE. So we show the process 
of the resonant transformation of the solution and the appearance of the wave 
packets with the different carriers. 

Earlier the resonant generation of periodic waves by a small external force 
was investigated by a computer simulation |5J. The phenomenon of the genera¬ 
tion and the scattering of the solitary waves in the case of nonlinear Schrodinger 
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equation was asymptotically investigated in 00 . The problem on the genera¬ 
tion of a solitary wave with a single carrier was solved in Here we use the 
proposed approach to study the multiphase case. 

The proposed approach is based on a local resonance phenomenon. The local 
resonance in linear ordinary differential equations was investigated in papers 
um- Later this phenomenon was investigated in partial differential equations 
for linear case in. and for weak nonlinear case mma 

The goal of this paper is the following: to demonstrate that the passage 
through the resonances allows to obtain the sequences of solitary packets of 
waves with the different carriers. The packets with different carriers do not 
interact. 

This paper has the following structure. The first section contains the state¬ 
ment of the problem, main result and example. The second section contains the 
asymptotic construction before the resonance. In the third section we construct 
an internal asymptotics in the neighborhood of the resonance curve. In the 
fourth section we construct the asymptotic expansion after the passage through 
the resonance. Five section contains the analysis of the resonances in high-order 
terms of the asymptotic solution. All asymptotics are matched. 

1 Main result 

Let us consider the Klein-Gordon equation with the cubic nonlinearity perturbed 
by external force with the finite collection of the different harmonics 

dfU — <9 2 U + U + yt/ 3 = e 2 ^ fk( ex ) ex P | — 1 + c.c., 0<£<1. 

fc=i l e J 

( 1 ) 

Here 7 = const; fk(y), k = 1,..., N and the phase function S(y, z) are smooth. 

We construct the formal asymptotic solution of the WKB-type. In the do¬ 
main when the first mode of the perturbation does not pass through the reso¬ 
nance the asymptotic solution has a form 



where 


lj(x 2 ,t 2 ) = j 2 (d t2 S) 2 - j 2 (d X2 S) 2 - 1, x 2 = e 2 x, t 2 = e 2 t. 

The leading-order term has an order of e 2 and oscillates. Such solution relates 
to the forced oscillations. 

On the curve lk{x 2 ,t 2 ) =0, k = N,..., 1 the frequencies of fc-th mode of 
the forced oscillations and a frequency of the eigen oscillations of the linearized 
Klein-Gordon equation are equal. Usually the curves with the such property are 
called resonant curves. The local resonance layer appears in the neighborhood 
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of the curve l k (x 2 ,h) = 0. Here we describe the passage through the resonant 
layer with the ordinary number k. 

After passage through the fc-th resonant layer a new eigen mode with the 
amplitude of the order e appears in the solution. It leads to changing of the 
number of harmonics in leading-order term. The amplitudes of the leading-order 
term satisfy (N — k + 1) nonlinear Schrodinger equations after passage through 
the fc-th resonance. The solution has a form 

N 1,0 

U = £^2 'i' $3 exp {i^j/e 2 } + 0(e 2 ), (3) 

j=k 

here the value k in © is defined by 


k = min j : lj > 0, j = 1,..., N 
and functions Ij(x 2 ,t 2 ) are arrayed: 

lj(x 2 ,t 2 ) < lm(x 2 ,t 2 ), l<j<m<N. 

The accurate formulation of the result for this paper is following 
Theorem 1 Let the asymptotic solution of m relates to forced oscillations: 

N f 

U = -e 2 p ex P {KjS)/e 2 } + 0(e 3 ) 

3=1 


in the domain In < 0 before the passage through the first resonance and the 
asymptotic solution has a form 

N bo 

U = e ^ \i> ^ exp{*<hj/e 2 } + 0(e 2 ), 

j=k+l 


in the domain —Ik £, lk+i £ before the k-th resonance curve Ik = 0. Here 

£ H b 10 = {=b$fc+l,±^ > fc+2, • • • ,i^Ar} 


1,0 


and the amplitudes satisfies (N — k) nonlinear Schrodinger equations 


1,0 


1,0 


1,0 


2 idt 2 ^jd t2 + d£. T ^ + *[<9 2 2 $j - <9 2 2 $,] 'f ^ + 

7 l ^ | 2 ^ = o, j = k + 1,..., N. 


(4) 


Then in the domain —Ik- 1 e, h £ the solution has the similar structure 

with the changes of the phase collection according to 

r% = r%u{±<s> k }. 

The number of equations in OF increases up to (N — k + 1). 
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Here we use the upper indexes b and p in our notations to specify the items 
in the domain before the passage through the k -th resonant curve and post 
resonance items. 

To illustrate the the main result let us consider equation 0 with the simplest 
driving force with two modes: 

F = /iOi) exp(*(t2/2 + x 2 )/e 2 ) + .M^i) expfy^ + 2x 2 )/e 2 ). 

In this case the curves of the local resonance are the lines t 2 = i/l7/16 and 
t 2 = y/2. In the domain t 2 < \J 17/16 the asymptotic solution is 

U = 0(e 2 ). 

The solution of this order with respect to e relates to forced oscillations. 

After the passage through the first resonance the solution is 

U = s vj/ $2 exp{*<f> 2 /e 2 } + 0(e 2 ), y/17/16 <t 2 <\[2 

1,0 

The amplitude T $ 2 satisfies the Cauchy problem for 

1,0 „ 1,0 1,0 1,0 

2 id t2 V $ 2 + df l6 V $ 2 +7| T $ 2 | : $ 2 = 0, 

variable = ti — \J\7/lQx\ , Initial condition is 

* $2l t2= yi77^ = /i(0(i + *)V^- 
After the passage through the second resonance the solution is 

U = e T exp{i$i/e } + e \j> $ 2 exp{i<f> 2 /e } + O(e ), t 2 > V2. 

1,0 1,0 

The amplitudes H/ and vj/ $ 2 are determined from the Cauchy problem for 
two nonlinear Shrodinger equation: 

1,0 „ 1,0 1,0 1,0 

2id t2 ^ +9| l4l ^ +7l =0, 

1,0 1,0 1,0 1,0 

2 id t2 $ 2 + <9f 26 ^ $ 2 + 7l T $ 2 |' f 3-2 = 0. 

Here £ 2 = t\ — \/2x\/2 1 initial conditions are 

1,0 1,0 1,0 

^ $2lt 2 =V2 = ^ ®2lt 2 = V2-0’ ^ 4 , ilt 2 = V2 = M&X 1 + *)V71. 

The solution of this Cauchy problem contains solitary waves if the initial data 
are sufficiently large m- 
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2 Pre-resonance expansion 

In this section we construct the formal asymptotic solution in the domain before 
the fc-th resonant curve. Here we use scaled variables Xj = e°x, tj = s :l t, j = 1, 2. 
The solution has a form 

OO 

U(x,t,e) = e* b 1 + Y,e n * b n > ( 5 ) 

n= 2 

N 10 

*? = E E exp{±*$ m (a:2,t2)/e 2 } T 

m=k -\-1 

The higher order terms are 

n—2 N , 

^n = E lni ( e ) E (E ^xp{dzZ^ > m (x2, ^ 2 ) I £ } \I/ ±«I> m 7 tl > ^ 2 ) “I - 

j=0 m—/c+1 ' ±<£ m 

X! ex p{*x(^ 2 ,i 2 )/e 2 } ¥ x (a;i,ti,t 2 )j, 

X^T bl nJ ' 

where • is a set of phase functions which is determined by 

T% = {±$ fc+1 , ±<& fc+2 , ■ • •, ±<M; ^ 2,0 = T? i0 U {±5,..., ±NS}, 

~^n,j = U Xm,ji + Xn 2 ,j 2 ~^Xn3,J3i Xn p ,j p € Y np j p . 

n\ + 712 + nz = n, 

jl + 32 + J3 = j 

r Y'bf 'V'b \ r v'b 

-*■ n,j n,j \ 1,0* 

This solution contains the two parts. The first part of the solution has the 
leading-order term of the order of e. The phase set of this part is Tj 0 = 
{±$fe+i, ±3>fc + 2, • • •, The modes have appeared in the solution due to 

the crossing of the previous (N — k) resonant curves lj =0. The terms of the 
order of e 2 with the phases ^ Tj 0 relate to the forced oscillations and describe 
the behaviour of the solution before the first resonance curve In = 0 where the 
forced oscillations take place only. 

The asymptotics constructed in this section is valid as lk+i e, — Ik £• 
The result of this section is formulated below. 

Note. Expansion © contains the terms e" ln m e. It looks a little bit un¬ 
expectedly because of the original equation does not contain the logarithmic 
terms with respect of e. But these terms naturally appear due to the passage 
through the previous resonant layers, see Lemma [S] in subsection Id. 3. II 

Let us construct the formal asymptotic solution for equation m in form ©. 
We substitute © in equation © and collect the terms of the same order of e. 
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As a result we obtain a recurrent sequence of equations for the coefficients of 
the asymptotics. 

Terms of the order of e give us the equations for the phase functions 

- (^= 2 ^j) 2 -1 = 0, j = k + 1,..., N. (6) 

Initial data for the phase functions <f>j, j = k + 2,... ,N are determined by their 
values on the curve lk+i = 0. Initial data for is related with 


^ik+i=o — (fc + l)'S'k+ 1 =0j 0t2$jk+i=o — (& + 1)^*2 ‘S'k+i=o- (7) 

The terms of the order of e 2 give us the homogeneous transport equation 


d t2 $jd t i T <r>, - d X3 ®jd Xl T = 0. j = k + l,...,N. (8) 

This equation allows us to determine the dependence of the leading-order term 
on characteristic variable £. Equation © along the characteristics 


dx i 

dQ 


- 9x2 5 


dti 

dQ 


= dt2 


(9) 


can be written in the form of ordinary differential equation 

1,0 


d v]/ $. 


( 10 ) 


1,0 


It yields $, depends on A, where the A is defined by 


dx\ u ^ ^ 

d & ' 2 '■ d & '" 2 '■ 


Among the terms of the order e 3 we collect the terms which oscillate as 
exp (i&j/e 2 ). It gives 


f 2,0 2,0 \ 

2 i ld t2 $jd tl ^ - d xa ®jd xi 4-^1 + 

2id ta <f> j d ta i ^ + [{d tl Q 2 - (d xi ^) 2 ]dl^ T ^ + 

’] ^ +7l ^ -fj 2 'h <s 3 - = 0. 


It is convenient to write this equation in the form of ordinary differential equa¬ 
tion in terms of characteristic variables 


2,0 

d T $, 

^cT 


1,0 


1,0 


= - 2 idt 2 $jd t2 # 9j - [(dt^j) - (d Xl &) ]<% c . $ 
„ „ 1,0 1,0 1,0 


( 11 ) 
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Equation El shows that the right hand side of equation ra does not depend 
on Q. To construct the bounded solution of GH> we demand the right hand 
side of the equation is equal to zero. It allows to determine the dependence of 
the leading-order term on slow variable t 2 ■ The amplitudes of the leading-order 
terms satisfy (N — k) nonlinear Schrodinger equations 

2 idt a $jdt 3 T + [(9 tl ^) 2 - <10 + 

+i[dt 2 ®j - i ^ +t| ^ $ 3 | 2 ^ 4> 3 - = 0, j = k + l,...,N (12) 


1 ’° 

The initial conditions for T 3 ^ are 


1,0 

^ < ^ > fc+ibfc+i =0 


1,0 1,0 

4' $jk=o = 4» 4>i|ifc=+ 0 ) j = k + 2,...,N] 


f 

J —( 


dcrfk+i{xi) exp(i 


d/-L\k+i(xi,ti,e)). 


Integration in this integral is realized in the line of characteristic direction con¬ 
nected with ©• 

The equations for the higher-order terms are obtained by the same manner 

/ n+l,fc n+l,fc \ n,k n,k 

2 ild t 2 $ j d tl 4- - d X 2 $jd Xl T j =2 id t 2 $jd t2 T 3., i>t Ai 41 


n,k n—l,k 

V Zj+dt^jd L T $,-7 


v ki,k2 h,h mi,m 2 

2 T a T 0 T 6 , 

k 1 ,k2,li,l2, mi,m 2 , cz,(3, <5 


where fci + h + TOi = n + 2, k 2 + h + rn 2 = k, a + (3 + 5 = 4>j, a £ /3 £ 

T b icT 6 

1 /i,Z 2 ’ a ^ i mi,m 2 ' 

To construct the uniform asymptotic expansion with respect to Q we obtain 
the linearized Schrodinger equation for higher-order term 

n,k n,k n,k 

2 idt 2 $jdt, 4- <i>, • il Ai 4- *, + i[d ta $j - d xa $j] 4- = 

n—l,k x ^ ki,k2 h,h mi,m2 

-<9 tl £df. t2 4-^-7 2^ 4- 4-/? 4- 5 , (13) 

ki,k2,li,l2 ,m i ,m2 ,a,( 3 ,S 


where k\ + 1\ + mi = n + 2, k 2 + h + m 2 = k, a + f3 + 5 = , a £ fc2 , [3 £ 

T b A <= T b 
1 Zi,Z 2 ’ 0 1 m\,m 2 ’ 

n,k . 

The amplitudes $ x as x G T are determined by algebraic equations 

r o on n ik n,k 

[-(Xt 2 ) +(Xx 2 ) + 1 ] 4> X =E X , x¥=±®j, j = l,...,k. (14) 

Here the right hand side of the equation depends on previous terms and their 
derivatives 

n,k n—1,/e n— 1,/e n—2,k n—2,k 

F x = ~^Xt 2 9t 1 4 1 x 3~2iXx 2 dxi 4 f x ~^Xt 2 ^t 2 4 f x ~^[Xt 2 t 2 ~ Xx 2 x 2 ] 4 f x ~ 


7 


n—3,k n—4,k 

dlt 2 * v -9, 


v fei,fe2 h,h mi,m2 

t 2 t 2 ^ x “7 2^ ’f a ^ 0 ^ 5, (15) 

ki,k2,li,l2 ,m 1 ,m2 ,a,f 3,8 


where fci + Zi + mi = n — 4, k 2 + Z 2 + m 2 = k, a + /3 + 5 = x, «£ T£ fca , /3 e 

T b A p t 6 
1 /i,Z 2 ’ a ^ 1 mi, m2' 

Note. It’s necessary to note that the multiplier [—(xt 2 ) 2 + (Xx 2 ) 2 + l] can 
be vanished on some curves. It leads to the resonances for the higher order 
terms of the asymptotics. The passage through the resonances does not change 
the leading-order terms of the asymptotic solution. We discuss this passage in 
section 0 

n,k 

In this section we pay a special attention on the amplitudes vj> x with the 
phase function kS. These amplitudes have the strongest order singularity on 
the curve l k = 0. The multiplier [— (xt 2 ) 2 + ( Xx 2 ) 2 + l] on the left hand of 
equation m equals zero on the curve Ik = 0. Here we explicitly write out a 
few first correction terms: 

fk 


2,0 

^ kS= -y-, 

Ik 

n ;9 xl f k d X2 S 

Vks= 2 ik - -2 -. 

l k 


(16) 

(17) 


4,0 2 ikf k [d t2 Sd t2 l k - d X2 Sd X2 l k ] - 4 k 2 (d X2 S) 2 d*J k 

4'fcS=-—---——--- ! - 

L k 

2ikd t2 f k d t2 S + dlj k + ikdf 2 Sf k 
l 2 

l k 


(18) 


The formula for the n-th order term has the form 


n,m 1 n—4,m n—2,m n—2,m 

C\Z 


T ks= 7 - 

f'k - 


d t2 T kS +2 ikd t2 Sd t2 T kS +ikS t2 t 2 T kS - 


n—2,m n—2,m 

2 ikd X2 Sd X2 % kS ~ikdl 0 S % k s ~df 


n—3,m 


^2 


n— 4,m 


n—2,m 

T kS - 

n—l,m 


2d 2 1X2 ks ~d X2 T ks -2ikd X2 Sd Xl T k s + 


7 


E 


ni,mi ri2,m2 n^,m^ 

T fci HI fc 2 HI fe 3 


ni+ri2+n3=n, 

ki+k2+k$=kS 

7711 +7712 + 777 - 3=771 


(19) 


The main order of the singularity corresponds to the terms when m = 0. 
This order depends on the number of correction term. 







n,m 

Lemma 1 The coefficient (f E T\ 0 U {kS} has the following behaviour 

T„= 0(Z- ( —- 1} ), l k - -0, (20) 

Proof. At first we prove formula for the phase <p = kS. The validity of 
formula in for n = 2,3,4 and m = 0 directly obtains from EO), EJ, m- 

n- 1,0 

Suppose now that this formula is valid for the term vj/ tp . The increase of the 
order of the singularity as l k —> 0 takes place due to differentiation with respect 
to X 2 , t 2 and the nonlinear term in formula m- Differentiation of the terms in 
formula U3 leads to formula m for m = 0 . 

The validity for values m > 1 and other values of <p can be obtained by 
direct calculations as it was shown above. Lemma is proved. 

These lemma allows us to write the asymptotic representation for the coef¬ 
ficients 

OO 

Tb 7TL 7b 7Tb 

4 „= Y, H l l V e T 1,0 u {fc 5 }, l k -> -0. (21) 

j=— (n—2m— 1) 


n,m 

The terms ^ 0 U{kS} have the more weaker order of the singularity 

as lk —> —0. Equations m for these values of tp have the same structure but 
the multiplier 7 - is regular as l k —* — 0 . 

L j 

The domain of the validity as l k —* —0 for the formal asymptotic solution in 
the form © follows from the relation 


e^n+l 

ijfb 


< 1. 


It yields 


The following theorem is proved. 


-Ik » £■ 


Theorem 2 In the domain —l k 3> £, l k +1 £ the formal asymptotic solution 

of equation m has form The coefficients of the asymptotics are defined 

either from differential equations mi, mi or from algebraic equations 

(7?|). nk>- 


3 Internal asymptotics 

This part of the paper contains the asymptotic construction of the solution for 
equation m in the neighborhood of the curve l k = 0. The domain of validity 
of this internal asymptotics intersects with domain of validity of expansion ©. 
These expansions are matched. Here we present asymptotic constructions for 
arbitrary k in details. It is necessary to note this section also describes the 
transition from the forced oscillations of 0(e 2 ) to the solution of 0(e) when 
k = N. 
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Let us construct the internal asymptotic expansion in the domain \lk\ <C 1. 
Denote 

A fc (xi,ti,£) =-/ fc (ea:i,eti)- (22) 

£ 

Theorem 3 In the domain \lk\ 1 the formal asymptotic solution for equation 
m has the form 


U(x,t,e) = £W 1 + '^£ n W n , (23) 

n =2 


where 


1,0 


N 


Wi =WkS exp{ikS(t 2 ,x 2 )/e 2 }+ E E exp{±i<& m (x 2 ,t 2 ) / £ 2 } 


m=k +1 ±$ t ? 


1,0 

W ±$ m (^l,tl,t2) 


W n = y ln J (e) WkS exp (ifcEiiiE.) + 
j=0 £ 

n —2 / N 

E lnJ ( £ )( E E ex P{±*$m(z 2 ,i 2 )/e 2 } W ±* m (xi,t 1 ,t 2 ) + 

j—0 ^ m—k-\-l ±«I> m 

2 n ’-? \ 

2^ exp{*x( ;r 2 ) t 2 )/£“} W x (xi,ti,t 2 ) J, 

xer’ nJ J 

where T' n ■ is the set of phase functions which is determined by 
Ti i0 = ±$ fc+ i. ±$ fc+2 ,..., ±$jv, ±kS ; 

Y 2 ,0 = ±$ fc+ i, ±$ fc+2 ,..., i'Ljv, ±5,..., ±NS , 

~^n,j = U Xmji + Xn 2 ,j 2 + Xn 3 ,j 3 i Xn p ,j p € T n p ,j p ■ 

ni +n 2 + ri 3 = n, 

ji + J2 + J3 = j 

T / r Y' \ 'Y' 

n,j — 1 n,j \ 1 1,0- 

n.j ■ 

TTie function WkS is solution of the problem for equation with zero condi¬ 
tion as Afc —> —oo. 

There is a difference between asymptotics (Pl> and external pre-resonance 
asymptotics 0. Asymptotic expansion defines the generation of the k- 
th mode by the local resonance for the leading-order term of the asymptotic 
expansion. The difference is the collection of phase function. The collection of 
phases Tpo contains (N — k + 1) phases but the phase ±<f>fc does not exist in 


10 



the k -th resonant layer. There phases ±kS are presented in the T^o and T 2 j o 
collections. 

The proof of theorem [3] consists in three steps. First we derive equations 
for the coefficients of the asymptotics. Second we solve the problems for the 
coefficients. And third we determine the domain of the validity for expansion 

ed- 

3.1 The equations for coefficients 

In the domain 1 <C Afc <C e _1 both asymptotics © and (El are valid. This 
fact allows us to obtain the asymptotic representation for coefficients of ED- 
Substitute h = sXk in formula ED and expand the obtained expression with 
respect to powers of e. It yields 

W v = E X~ k jn %’ m Ux 2 ,h,x 1 ), Afc —> —oo. (24) 

j=n—2m— 1 

Let us obtain the differential equations for the coefficients of asymptotics 
ED- Substitute ED in equation © and collect the terms with equal powers 

n,j 

of small parameter and exponents. It yields the equations for coefficients W<p, 
tp £ Y raj -. In particularly, the terms of the order e 2 give us the equations for the 
amplitudes of the leading-order terms 

2ikd t2 Sd tl Wks -2 ikd X2 Sd Xl Wks -A fc WkS= fk, (25) 


1,0 


1,0 


1,0 


2 id t2 $jd tl W<s>j -2 id X2 $jd xl W<s> 2 ~h j (x2,t 2 ) W<s>j= 0, j = k + 1 ,..., N. 


(26) 


1,0 


and complex conjugated equation for W-<p, P £ Tfi,o(^)- 

The relation of the order e 3 in equation © gives 2 (N — k) + 2 equations: 


2,0 


2,0 


2,0 


1,0 


1,0 


2 ikd t2 Sd tl Wks -2 ikd X2 Sd Xl Wks -A k Wks= d Xl Wks -df, Wks ~ 


1,0 


1,0 


1,0 


-ik[df 2 S - d z X2 S} Wks -2 ikd t2 Sd t2 Wks +2 ikd X2 Sd X2 Wks 


1,0 


N 

-3 7 \\WksY+ E I 

m=k -\-1 


1,0 


1,0 

Wks ■ 


(27) 


2,0 


2,0 


2,0 


2 id t2 $jd tl W<f>j -2 id X2 $jd Xl W~h j (x 2 ,t 2 ) W$ j = 


1,0 


1,0 


1,0 


d Xl W^j -d$ ± W^ j - -dl 2 $j] W&j -2idt 2 $jd t2 W^ j +2 id X2 Sd X2 


-37 


1,0 

WkS 


N 

E 

m=k -\-1 


1,0 


1,0 


W ■ (28) 
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The higher-order terms are calculated by the same way. In particularly 

n,j 

W v ,p £ lyo is determined by differential equations. Here we represent the 
equation for amplitude under = kS. The equations for other phases ip £ Tqo 
can be obtained in a similar. 


n,j n,j n,j n,j 

2 ikdt a Sd tl Wks -2ikd X2 Sd Xl Wks -A k Wks=F k ■ 


(29) 


The right hand side of equation has the form 

n,j n—l,j n—l,j 


n- l,j 


F k = -2ikd t2 Sd t2 W kS +2ikd X2 Sd X2 W k s +k{d t2 S) 2 W kS - 


o 0 n— l,j „ n—l,j _ n—l,j n—2,j n—2,j 

k 2 {d X2 S) 2 W ks -dl W ks +d 2 Xl W ks ~d t2 d tl W ks +d X2 d Xl W ks - 


n-3,j n-3 ,j 

d 2 2 W ks +C W ks -7 


E 

ni + n 2 + ri 3 = n + 1, 
ji + 3 2 + J3 = 3, 
X1+X2+ X3 = kS, 
Xj £ F nj {k), j = 1,2,3 


nidi n 2 j 2 n 3 ,j 3 

W X1 W X2 W X3 . (30) 


n,k 

The amplitudes W x as x Yi.o are determined by algebraic equations 


,k n.k 


[ _ (x*2) 2 + {Xx 2 ) 2 + l] W x — F' 


(31) 


Here the right hand side of the equation depends on previous terms and their 
derivatives 


n,k n—l,k n—l,k n—2,k n—2,k 

F X = ^Xt2^ti W X ~X2iXx 2^X1 W X — ^X*2^*2 XV X [Xt2*2 _ XX2X2 ] XV X — 


dlt 2 


n—3,k 

W y -d 2 


*2*2 


n—4,/c 

W x -7 


E 


N,k 2 ,h,h j TTb i ,m 2 ,a,/3,6 


N , k 2 mi,7712 

WaW /3 W 6 , 


(32) 


where N + l\ + mi = n — 4, k 2 + l 2 + m 2 = k, a + f3 + 5 = x, a £ T nm 1 P £ 

"£i lt i 2 , S £ T mijm2 . 


3.2 The solvability of equations for higher-order terms 

n,j 

In this section we present the explicit form for the higher-order term Wks and 
investigate the asymptotic behaviour as \ k —> ± 00 . 
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3.2.1 Characteristic variables 


n,j 

The function WkS satisfies equation The solution is constructed by char¬ 
acteristic method. Define the characteristic variables a, £. We choose a point 
(xi,t°) such that d X2 lk\( x ojl) ^ 0 as origin and denote by a the variable along 
the characteristic family for equation We suppose a = 0 on the curve 

Afc = 0. The variable £ mensurates the distance along the curve Afc = 0 from 
the point This point (x®, t ( {) corresponds to £ = 0. Let the positive 

direction for parameter £ coincides with the positive direction of x -2 in the neigh¬ 
borhood of (x°, ti). 

The characteristic equations for (l29l) have a form 

= 2kdt 2 S(sxi,£ti), Pp- = —2kd X2 S(sxi,eti). (33) 

da da 

The initial conditions for the equations are 

£i|<t=o = 2u, ii| CT =o=i?- (34) 


Lemma 2 The Cauchy problem El. El for characteristics has a solutions 
as \a\ < Ci£ _1 , ci = const > 0. 

Proof. The Cauchy problem (1551) . (1551) is equivalent to the system of the integral 
equations 


— 2 f kd X2 S(sxi,eti)d(. (35) 

J o 

Substitute £2 = (ti — ti)e, X 2 = (aq — x\)e. It yields 

pea pea 

t 2 = 2 / kd t 2 S(x2—£X ( l,t2—x 2 = -2 / kd X2 S{x2—£x\,t2-et\)d^. 

Jo Jo 

The integrands are smooth and bounded functions on the plane Xz,t 2 ■ There 
exists the constant ci = const > 0 such that the integral operator is the con¬ 
traction operator as e|er| < ci. Lennna|5]is proved. 

Assumption. We assume that the change of variables (xi,ti) —> (<r, £) is 
unique in the neighborhood of the curve Afc = 0. This assumption means that 
the characteristics for equation lp|) do not touch the curve Afc = 0. It means 


ti = t\ + 2 f kd t2 S{ex\,et{)dC,, x\ = x® 
Jo 


d X2 l k d X2 S - d t2 l k d t2 S ± 0. 


It is convenient to use the following asymptotic formulas for change of vari¬ 
ables (xi,ti) —> (cr,£). 
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Lemma 3 In the domain |a| <C £ 1 the asymptotics as e —> 0 of the solutions 
for Cauchy problem Ei, O /lave the form 

N 

Xi(a, £, e) — Xi(£) = —2crkd X2 S + 2 e"a" +1 (? n (ea:i, eii) + 0(e N+1 a N+2 ), (36) 

n=l 

N 

t! (ff,£,e) -i?(0 = 2(7jfca t3 S + 2^e n <r n+1 /i ri (ea: 1 ,et 1 ) + 0(e JV+1 <7 jV+2 ), (37) 

n=l 


where 


9n — & ^ n (^X2 


cr=0 


1 K = k d^^ s) 


CT —0 


The lemma proves by integration by parts of equations 
The next proposition gives us the asymptotic formula which connects vari¬ 
ables a and A*, as a, X k —> ±oo. 


Lemma 4 Let be a -C e 1 , then: 


Afc = (p(f)<J+ 0{ea 2 ), $(£) = 


dX k 

da 


a —> oo. 


a—0 


Proof. From formula we obtain the representation in the form 

OO 

Afc = X! A£(aq^i^cr-V -1 , 

3 =1 


where 


It yields 


Let be 


1 d? 

Ai(xi,fi,e) = -t — A fe (a:i,ti,£)| <T= o- 
j! acr-? 

\ ^A fe I , 2 ^ Afc 

Afe = -7— I _ n cr + 0(ea -- 77 - 
da lCT -° v da 2 

d 2 lk 


)• 


da 2 


> const, f £ R. 


The function dA& /dcr is not equal to zero 

^ = \( ~ kd X2 X k d X2 S + kd t2 X k d t2 S ) ^ 0. 


dX k 


da 2 

Let us suppose dX k /da > 0. It yields 


Afc = (j)(£)<J + 0 (£<j 2 ), = 


da 


CT = 0 


The lemma is proved. 
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3.2.2 Solutions of the equations for higher-order terms 

n,j 

The higher-order terms W±kS are solutions of equation m with the given 
asymptotic behaviour Afc —> — oo. Equation can be written in terms of 
characteristic variables as 


d n d n,j n,j 

i~r- WkS ~Xk WkS=FkS ■ 
da 


(38) 


Lemma 5 The solution of equation VM) with the asymptotic behaviour as 
Xk —► —oo has a form 


n,j 


n,3 


Wks=exp(-i d(\k{xi,ti,e)) / d(Fks(x 1 ,t 1 ,£)exp(-i dxXk{xi,ti,e)). 
J 0 J — oo J 0 

(39) 

Proof. By direct substitution we see that expression (1351) is the solution 
of l!331) . The asymptotics of this solution as Afc —> —oo can be obtained by 
integration by parts and substitution 


da 


= 2 kd t2 Sd tl - 2 kd X2 Sd Xl . 


It yields 


n i3 v—"v 

Wks= y , 

j= o 


2 kd t2 sd tl - 2 kd X2 sd Xl ' 

iXk 



_ rl i3 

Fk 

) 

iXk 


, X k — i ► -oo. (40) 


From formula CT we obtain that formulas m and (El are equivalent. The 
lemma is proved. 


3.3 Asymptotics as —> oo and domain of validity of the 

internal asymptotics 

The domain of validity of the internal expansion is determined by the asymp¬ 
totics of higher-order terms. In this section we show that the n— th order term 
of the asymptotic solution increases as A)! -1 when Afc —> oo. This increase of 
higher-order terms allows us to determine the domain of validity for internal 
asymptotics El as Afc —> oo. 

3.3.1 Asymptotics of higher-order terms 

This section contains two propositions concerning asymptotic behaviour as A k ~► 
oo for higher-order terms in El- The first lemma describes the asymptotic 
behaviour of higher-order terms as Afc —> oo and the second one contains a 
result about asymptotics of the phase function. 
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Lemma 6 


n,j 

The asymptotic behaviour of WkS as 1 X k £“ 


has a form 


n ' 71-1 i~i 

WkS = EE( A*. In” 1 1Afe| W'fe’ m) (0 ) exp(—z j d(Xk(xi,ti,e)) + 

j—0 m =0 


+ E 


2 kd t2 Sd tl - 2 kd X2 Sd xl ' 
iX k 


m 

m. 7 

' FkS 

) 

iXk 


m =0 

Proof. Let us calculate the asymptotics of the leading-order term 

d(fk( x i) exp(i 

to 


(41) 


10 r<r /< !■<? 

Wks= exp(-i / dCA fe (a;i, ti, e)) / dCf k (xi)exp(i dxXk(xi,ti,e)) = 
J 0 J—co J 0 


exp(— i 

/ d(X k {xi,h,£)) I 

nOO 

' d(f k (x i)exp(z 

[ dxX k (xi,ti,£))- 

j 

'o J 

— OO 

1 0 

exp(— i 

/ d(X k (xi,ti,e)) 
Jo 

/‘CO 

j dC,fk{xi) exp(« 

[ dxX k (xi,ti,s)). 
Jo 


Further by integration by parts of the last term we obtain formula SB as n = 1, 
where 

i o r°° 

Wks\0= daf k (xi)exp(i / dxA fe (a;i,ii,e)), 

J—co J 0 

1,0 

F ks— fk{x i). 

2,0 

To calculate the asymptotics of WkS in formula we use the asymptotics 
with respect to a of the leading-order term. Integral ilMl) contains the term 
with linear increase with respect to a when n = 2. We eliminate this growing 

part from integral explicitly. The residuary integral converges as <r —> oo. It 

1,0 

can be calculated in the same manner as it was calculated for WkS- It yields 
formula m as n = 2, where 


2,0 

W 


(1,0) 

kS 


1,0 

(0 =W 


(0,0) 

kS 


( 0 - 


The same direct calculations are realized for the n—th order term. The lemma 
is proved. 

To complete the proof of theorem [3] we need to obtain the domain of validity 
of asymptotics G3 • The formal series 423 is asymptotic when 


eW n+ 1 
W n 


< 1 , 


0 . 


Lemma ED gives Afc -C s 1 . After substitution X k = eh we obtain (j C 1. 
Theorem |3 is proved. 
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3.3.2 Asymptotics of the phase function as Afe —» oo 

To obtain the asymptotics as A^ —> oo we need to derive the asymptotics of the 
phase function in formula 6D- 


Lemma 7 ^4s X k —> oo: 

[ d£\k =-^2 + -{d Xa S(xi-x°) + dt a S(ti-t°))+0(e\l). (42) 

Jo £ £ 

Proof. Substitute the asymptotics of A k from lemma ED Calculate the 
asymptotics of the integral in formula 03 


dCXk(xi,ti,e) 



(-kd x J k d X2 S + kd t2 l k d t2 S) C + 0(£ c 2 ) 


(- kd X2 l k d X2 S + kd t2 l k dt 2 S)£- + 0(£ct 3 ). 

The asymptotics of the phase function kS(x 2 ,t 2 ) in the neighborhood of the 
curve l\ = 0 is represented by a segment of the Taylor series. It yields 

UQ | 

~2 = -(kd X2 S(x! - x°) + kd t2 S(t\ - t°))+ 

£ z £ 

— {kS X2X2 (xi — .t^) 2 + 2 kS X2 t 2 ( X\ — Xi)(ti — 1 5) + kSt 2 t 2 (ii — ^i) 2 ) + 

o(£(\t 1 -t° 1 \+\t 1 -t° 1 \n 

Substitute instead of (x\ — X®) and (t\ — tf) their asymptotic behaviour with 
respect to £ from lemma |3] This substitution and result of lemma 0] complete 
the proof of lemma [7| 

The asymptotics as At —> — oo contains fast oscillating terms with phase 
functions <hfc + i,..., mS,m £ Z, m ^ k. The leading-order term of the 
asymptotics as X k —> oo contains the oscillations with an additional phase func¬ 
tion. We obtain this result from lemma El Denote this new phase function by 
•Pfc(a^ 2 , h)/£ 2 - The asymptotics of this function is obtained in lemma 0 The 
nonlinearity and additional phase function lead to more complicated structure 
of the phase set for higher-order terms of the asymptotics as X k —> oo. 

Lemma 8 The phase set T n j for the n—th order term of the asymptotics as 
Afc —> oo is determined hy formula 

Tyo = {±$fc, • ■ ■ , ±4>Ar}; T 2 ,o = {±$fc, .... id* at, ±S ,..., ±NS}, 

*-n,j = U.Xniji + Xn 2 ,j 2 + Xn 3t j 3 i Xn m ,jm ^ 
where n\ + n 2 + n 3 = n, j\ + j 2 + jo = j- 
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The proof of this lemma follows from the asymptotic formula for n —th order 
term. Representation formula GH> and lemma [B] allow us to construct the 
asymptotics as A —> oo of the internal expansion in an explicit form 


exp 


N n— 2 / n— 1 n— 2 

17 = E £ ” E lnP ( £ ) lnm i a i w^ m \o 

n= 1 p—0 ' j— 0 m= 0 

.(k 


- i[ ~{d X2 S{x i - x") + d t2 S(ti - tp) + 0(sA J ) 


+ 


N n —2 / oo 

n —1 p =0 ^ j =0 


2fe(9t 2 59t! - 2kd X2 Sd Xl 
iX 


n iP 

F ks 


iX 


exp 




iV n—2 ✓ 

E] e ” El lnP ( £ ) ( EE ex p 

«=2 p=0 A v6Tf„, p \feS 


f.(p(t 2 ,a; 2 ) 

nr -> j- - 


c.c. (43) 


This representation and formula 1321 complete the proof of the lemma. 


4 Post-resonance expansion 

In this section we show how to connect the external asymptotic solution after 
the passage through the resonance with the solution in the neighborhood of the 
curve Ik = 0. 

There is a difference between these solution. The leading-order term of the 
solution in this section includes an additional mode. This new mode relates to 
the phase which was generated by the passage through the resonance near 
the curve Ik = 0. It gives 


T? i0 = T?, 0 U {±d>,} (44) 

This section consists in two parts. The first part contains the construc¬ 
tion of the formal asymptotic solution. This solution is similar to the solution 
from section |2] Asymptotic behaviour for higher-order terms as Ik —> 0 follows 
from section rnoi In the second part of this section we determine the domain 
of validity for this external asymptotics near resonance curve Zfc(a; 2 ,f 2 ) = 0. 
The matching method gives us the initial conditions for the coefficients of the 
asymptotics. 

4.1 Structure of the second external asymptotics 

Let us construct the formal asymptotic solution of the following form with the 
changing of the phase set according to (TH1) . 


U(x,t,e) = e4'(’ + <45) 

n=2 
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Y E exp{±i$ m (x 2 ,t 2 )/e 2 } <5 ±& m {xi,t 1 ,t 2 ) 

m=k +1 rt^m 


IL - / 

^ n E ln J (e) E E exp{±i$ m (a; 2 ,t 2 )/e 2 } 'S' ±<s> m (xi,h,t 2 ) + 

j =0 m=k +1 ' ±&m 

X—^ 9 n ’J \ 

E exp{ix(® 2 ,i 2 )/e } 'h x {x 1 ,t 1 ,t 2 ) I, 

X€l p ' n ,i 

where T p, nj is a set of phase functions which is determined by 

T? |0 = {±$ fc , ±$fc+i,..., ±<M; Yo = T lo U {±5, ■ • •, ±NS}, 

^n,j = U Xni,ji +Xn 2 ,J 2 + Xn 3 ,j 3 j Xn p ,j p € 

ni + ri 2 + 7 x 3 = n, 
ji + h + J3 = j 

T p, nJ = T>j\{±*k, ±*k+i, • • • - ±<M- 

Substitute this asymptotic solution in original equation and collect the terms 
of the same order with respect to s. After collecting the terms with the same 
phase functions we obtain the recurrent system of equations for the coefficients. 

The terms of the order e 1 give us the equation © for the phase function 
of eigen oscillations. The initial data is determined by matching condition 0 
with the changing of indexes from (fc + 1) to k. 

The difference between this asymptotic constructions and constructions from 
section |3 consists in changing of the number of terms of the order of e. 

1 ’° 

The amplitudes vf $ are determined from 

2 idtz&jdt, T <iq + 

+i[dt 2 ®j - d 2 2 $j] i +7l ^ *il' S ^ = 0, j = k,...,N (46) 

The higher order terms of m satisfy equation O- Here we present the 
following lemma which allows to match solution 63 and solution iTziil) . 

Lemma 9 The asymptotics as Ik —> 0 of the solution of equation m has the 
form 


n,p 

'h <s> k (xi,h,t 2 ) 


i -u-i) n 

E E lU^r+od), h-> 0 - 

j=— (n—2p— 2) m=0 

(47) 
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Proof. Determine the order of the singularity of the right hand side of the 

equation as l k —> 0. First consider equation m for n = 3,p = 0. The solution 

3,0 

of this equation gives us the coefficient vj/ . The nonlinearity contains the 

2,0 1,0 2,0 

term | k s | 2 . The function has the singularity of the order l k 

as lk —■> 0. It determines the order of singularity for right hand side l k 2 . We 

3,0 

construct the asymptotics of T x in the form 


3,0 3,0 

'F = >F 


— 1,0 7—1 
l k 


3.0 

+ 'F 


0,1 


3,0 

ln(4)+ 'F 


1,1 


3,0 

lk ln(4) + 


(48) 


Substitute (J3HJ in equation for n = 3. It leads to recurrent system of equations 
3,0 / • \ 

for coefficients <F 


3,o 1,0 2,0 „ „ 

-2id t2 <f> k dt 2 lk 4- L ’ = - 'F s\ l k , 


'S’k 

0 £ 1J = i[ 3 i° i-«], 

2 id t2 $ k d t2 l k i ^ = L[<F i 0 ^]. 

Here we denote the linear operator by 

£[*] = 2id t2 $kdt^ + dl^+i[dl$>k-dl 2 $k]*+l(% | 2 4' + ( 1 i°$J 2 vF*). 

3(0 

The regular part \F of the asymptotics satisfies the nonhomogeneous 
linear Schrodinger equation. The right hand side of the equation is smooth 


3,0 

F / ['F$ fc ] 


3,0 /1 1 \ 

~h In |Zfc|T['F & 1 '] 


2 id t2 $kd t2 l k 


3,0 

>F 


( 1 , 1 ) 


The initial condition for the regular part of the asymptotics is determined below 
by matching with the internal asymptotic expansion. 

n ,p 

The structure of the terms \F for n > 3 has a similar form. The right 
hand side of equation m depends on junior terms. These singularities can be 
eliminate 


n,p 

F $ k = 


— (n—2) -j-\-1 

E E 




n,p 

f k 


(j,m) 

$1 


n,p 

+ F k $ 1 . 


j = 0 m=0 


n ,P , ■ , 

The coefficients f k do not contain singularities as l k —> 0. These coefficients 
are easy calculated. 

The direct substitution of <G3 in equation and collecting the terms with the 
same order of l k complete the proof of lemma El 
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4.2 The domain of validity of the second external asymp¬ 
totics as Ik —■► +0 and matching procedure 

The domain of validity of the second external asymptotics is determined by 


zK+i 

*5 


< i. 


Formulas (TTTjll and m give the condition 


Ik » e. 


The domain \lk <C 1 of validity of the internal asymptotics and domain of 
validity of the second external asymptotics are intersected. This fact allows 
to complete the construction of the second external asymptotics by matching 
method ca The structure of singular parts of the internal asymptotics as 
Afc —> Too and external asymptotics as Ik —> 0 are equivalent. The coefficients 
are coincided due to our constructions. The matching of regular parts of these 
asymptotics takes place due to 


n.i 


k=o =w ks 0) (0- 


The function W ^s°\0 is determined in lemma HO 


1,0 


In particular, the initial condition for the amplitude <$ fc 
the new generated phase has a form 


1,0 

T 4>fek=0 = 


/ OO n<7 

dcrfk(xi) exp(* / dxXk(xi,h, 

-oo J 0 


corresponding to 

e)). 


1,0 

The initial data for the others amplitudes vj/ j = k + oi the 

leading-order term of (TTbli are represented by values of these amplitudes on the 
curve h = 0. 

The soliton theory for nonlinear Schrodinger equation leads us to the fact 
qo 

that the function <F ^ contains the solitary waves when is sufficiently 

large. 


5 Resonances in the higher order terms of the 
asymptotics. 

In this section we discuss the passage through the resonances of higher order 
terms of the asymptotic solution only. Such passage does not change the leading- 
order term and leads to changing of the solution in higher order terms only. 
Equations m are solvable while the multiplier 

I[x] = [ _ (xt 2 ) 2 + (xx 2 ) 2 +1] ¥= o. 
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Let us denote a particular value of x such that l[x] = 0 by Xr- In the neigh¬ 
borhood of the curve l[Xr] = 0, Xr G T Hy k the coefficients of the asymptotic 
expansion has a form 


n,k 

'I'Xr 


= 0(l~ 1 [ Xr ]), l[ X r]-> 0, 


(49) 


n+ ^\=0{l-^ +1 \xr ]), l[Xr}^ 0 (50) 

These formulas give the domain of validity of 

l[Xr) > £• (51) 

In the neighborhood of the curve l[Xr\ = 0 we use a new scaled variable \ Xr = 
l X r/ £ - The formal asymptotic solution is constructed in the form 

U(x,t,e) = U n - 2 (x,t,e) + U res (x,t,e). (52) 


The asymptotic solution has two parts. The first one does not depend on the 
scaled variable A Xr . But the second part of the solution depends on A Xr , because 
the resonance on the curve A Xr = 0 appears only in higher order terms of the 
asymptotics. The terms U n - 2 {x,t,e) and U res (x,t,£) have a form 


n— 2 m —2 / ^ 

exp{±*4>(a;2, t 2 )/£ 2 } 4> ±<s>(xi,h,t 2 ) + 

m —1 k —0 ' ±$ 

^^ m,k \ 

2 , exp{iXr(x 2 ,t 2 )/£‘ 2 } Xr (xi,ti,t 2 ) j. (53) 

cT' ' 


xrer' 


oo m —2 / rn k 

U res (x,t,e)= ^2 ^ ±4>(A Xr ,a;i,ti,t2) + 

m =n — 1 Zc—fl ' -J- cj) 

^ ^ 2 5 k "\ 

2^ exp{*Xr(® 2 ,t 2 )/£“} 4' Xr (A Xr ,a;i,ti,t 2 ) )• (54) 


m=n— 1 /c—0 


xreT; 


Substitution of (El into © gives the recurrent sequence of the problems for the 

m,k m,k 

coefficients of 15211 . Note that the coefficients 4 / ± $ and \J> Xr for m < (n— 1) 
are determined from the standard problems as shown in sectional 

m,k m,k 

The coefficients ±$ and T Xr for to > (n — 1) are determined as well as 
the internal expansion from section H Here we present the standard problems 
for the coefficients of the asymptotics without detailed derivation. The higher- 
order terms are calculated by the same way. In particularly, the terms in the 
case lower index is equal to Xr are determined by differential equations. 


m,k m,k m,k m,k 

2idt 2 Xrdt 1 4' Xr —2 id X2 Xrd Xl 4 1 Xr — A Xr 4) Xr = F Xr ■ 


(55) 
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m,k ___ ___ 

The right hand side F Xr of equation (1531) has form m with the changing of 
the index from x to \r- 

The terms in the case the lower index is not equal to Xr are determined by 
algebraic equations ED- 

The analysis of equation (1551) is realized in section EP1 The passage through 
this resonance layer leads to a new phase function appearance as was shown in 
We denote this new phase function by ip r . 

In the domain after passage through resonance layer the amplitude under 
this phase function is determined by 


n— l,k 


n—l,k 


2id t2 ip r dt 2 T Vr + ~ (d Xl Cr) ]9| r?r ^ Vr + 


n— l,k 1,0 _ n—l,k 

+i[d; 2 ifi r - d; 2 (fi r ] 4- v,r +71^4-1 I ^ 


= 0 , 


(56) 


here the variable is determined by 


dx\ dti 

d& = t2Xr ’ d^r 


9x 2 Xr ■ 


After the passage of the resonance layer the set 
following rule 


T 

77.— 


1 ,k 


1 — 1 ,k Xr 


Tn—i,fc is changed by the 


It leads to changing of the set of phases for higher order correction terms by the 
ordinary way. 
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